Angle-resolved photoemission studies of lattice polaron formation in the cuprate 
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To elucidate the nature of the single-particle excitations in the undoped parent cuprates, we have 
performed a detailed study of Ca2Cu02Cl2 using photoemission spectroscopy. The photoemission 
lineshapes of the lower Hubbard band are found to be well-described by a polaron model. By 
comparing the lineshape and temperature dependence of the lower Hubbard band with additional 
O 2p and Ca 3p states, we conclude that the dominant broadening mechanism arises from the 
interaction between the photohole and the lattice. The strength of this interaction was observed to 
be strongly anisotropic and may have important implications for the momentum dependence of the 
first doped hole states. 
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One of the most powerful aspects of angle-resolved 
photoemission spectroscopy (ARPES) is the potential 
to quantitatively compare experimental lineshapes with 
theoretical predictions, providing great insight into the 
nature of the single particle excitations [1]. For instance, 
good agreement between ARPES and theory has been 
demonstrated on model systems such as noble metal sur- 
face states 0, two-dimensional Fermi liquids H, and 
systems with strong electron- phonon coupling Ex- 
tending this quantitative analysis to the high-Tc cuprate 
superconductors has yielded a wealth of information re- 
garding the salient many-body interactions [J In 
principle, ARPES studies of the undoped Mott insulator 
should provide the ideal starting point to study the inter- 
actions in the cuprates, yet paradoxically, a logically con- 
sistent understanding of the spectra from the insulator 
has proven far more elusive than the doped compounds. 
However, an understanding of this simplest case is clearly 
essential before one can address the doped cuprates in a 
systematic manner. 

Although calculations based around the t — J model 
have successfully approximated the peak dispersion in 
the parent insulators @, Q, they have been unable to 
explain the broad lineshapes. Attempts to model spec- 
tra from the insulator within the same weak-coupling 
framework used for the optimally doped cuprates have 
proven unsuccessful. Very recently, a phenomenological 
picture based on Franck-Condon broadening (FCB) was 
proposed to explain these features. In this scenario, the 
broad peak was not a coherent quasiparticle (QP), but 
a manifold of states with the photohole strongly cou- 
pled to bosonic shake-off excitations. Within this model, 
previously unresolved issues such as the unusually broad 



width, the lineshape, and the large separation of the peak 
from the chemical potential, fi, have been explained in a 
self-consistent manner Q. However, the microscopic ori- 
gin of the broadening itself remains a fundamental but 
unresolved issue. In the past, most theoretical works 
considered only the effects of electronic correlations (i.e. 
Hubbard or t — J models). Recent work has indicated 
that electron-phonon interactions may also be significant, 
and thus it is important to determine the contribution of 
lattice interactions. Through a comparison with conven- 
tional "benchmark" states (O and Ca ip orbitals), 
with the complex and strongly correlated states of the 
lower Hubbard band (LHB), we conclude that lattice ef- 
fects play a substantial role in the observed broadening, 
implying that the photoholes form localized small po- 
laron states. In addition, we have observed a large mo- 
mentum dependence to the linewidth which may have im- 
plications on the anisotropy of the first doped QP states 
in the lightly doped cuprates. 

ARPES measurements were performed with a Scienta 
SES-200 electron analyzer at Beamline 5-4 of the Stan- 
ford Synchrotron Radiation Laboratory, as well as with 
a Scienta SES-2002 analyzer and a He plasma discharge 
lamp. Data were collected with photon energies between 
14-32 eV at the synchrotron, and both the He I (21.22 eV) 
and He II (40.8 eV) lines of the discharge lamp. Typical 
energy and angular resolutions were better than 15 meV 
and 0.35°, respectively. Multiple batches of Ca2Cu02Cl2 
were prepared using both a high pressure flux method Q 
and an ambient pressure technique (To| . 

In Fig. [T^, we show energy distribution curves (EDCs) 
from Ca2Cu02Cl2 taken near k = (7r/2, 7r/2) showing the 
LHB excitations. The chemical potential, /i, was deter- 
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FIG. 1: (a) The LHB and (b) O 2p.^ bands along (0, 0) - (vr, tt) 
taken with hv — 25.5 eV and T — 300 K, together with Gaus- 
sian fits (red) and orbital configurations (top). Lorentzian 
fits to the peaks at (7r/2,7r/2) and (tTjTt) are shown in dashed 
green, (c) Schematic of the Franck-Condon process, with the 
transitions from the initial state |^o^) to possible final states 
|*I'm~^) resulting in the observed photoemission intensity. 



mined from a Au reference in electrical contact with the 
sample, was found to be pinned within the charge- 
transfer gap by surface defects or impurities and varied 
from cleave to cleave. From a large series of measure- 
ments, we found that ^ was always separated from the 
LHB peak by at least 0.42 eV, and here we reference fi 
using this minimum value. Even at the top of the va- 
lence band, Flhb is extremely broad, while simple elec- 
tron scattering phase space constraints would predict ex- 
tremely sharp and long lived QP excitations. In addi- 
tion to the width, we also show in Fig. [1^ that the raw 
data can be fitted extremely well with a pure Gaussian 
function (with no additional parameters or background 
subtraction) over a wide range in both momentum and 
energy. At (7r/2, 7r/2), we also include a fit of the data 
to a pure Lorentzian lineshape (green dashed line) to il- 
lustrate the clear discrepancy between the experimen- 
tal data and the Lorentzian fit (which should closely ap- 
proximate the expectations from Fermi liquid theory near 
Fermi energy or at the top of the valence band). While 
the Gaussian lineshape is a much better fit to the data, it 
is inconsistent with a simple coherent single-particle spec- 
tral function, ^(k, w) On the other hand, a Gaussian 
lineshape is naturally expected for the simplest case of 
FCB in the limit of strong coupling, where the intensity 
is known to follow a Poisson distribution ^llj . 

In Fig. [TJ;, we show a schematic of the FCB process, 
where the initial state has finite transition probabilities 
to many closely spaced final states, and only the low- 
est transition would correspond to the QP. In principle, 
coupling of the electrons to any bosonic field could cause 
FCB, so studying the LHB alone is insufficient to dis- 



criminate between phonons or magnons. Therefore, we 
have compared the LHB with other states which are non- 
interacting with the spin system, but which necessarily 
inhabit the same lattice. We compare the LHB with 
the O 2p7r band, shown in Fig. [T}d, as first proposed 
by Pothuizen et al. in Ref. At k = (tt, tt) the O 2p.p; 
band is non-bonding with any Cu 3d states or the Zhang- 
Rice singlet state 12, 13] and thus is decoupled from the 
spin dynamics. For both the LHB and O 2p7r, we fit only 
the low binding energy half of the peak. This is due to 
the fact that both the LHB and O 2p^ states have neigh- 
boring states at higher binding ener gies , as shown by the 
cluster calculations of Eskes et al. [ij] and band struc- 
ture calculations [l^. Therefore, the spectral weight at 
higher energies (< -0.7 eV for the LHB and < -2.5 eV 
for the O 2pt^ band) contains contributions from multiple 
bands and therefore cannot be easily fit to a single func- 
tion. In our analysis, we are primarily concerned with 
the rapid suppression of the tail of spectral weight when 
going towards lower energies (closer to Ep), a hallmark 
of FCB, and which should not be strongly affected by 
higher energy states or background. 

The good agreement of the Gaussian fits with not only 
the LHB, but also the O 2p7r states indicates that a FCB 
lineshape is generic to all electronic states in the crys- 
tal. This would suggest that the photoholes form local- 
ized lattice polarons, since the O 2p^ states at (tt, tt) do 
not experience any coupling to the low-energy spin de- 
grees of freedom. A similar polaronic phenomenology has 
also been employed to explain the photoemission spectra 
from systems ranging from ionic insulators [T5|. charge 
density wave materials (K0.3M0O3, (TaSe4)2l) UJ, lUl, 
the Mott-Hubbard system VO2 [18], and the mangan- 
ites [1^, [2^ |2l| . The importance of lattice trapping in 
the Mott insulator becomes clear when considering that 
photoemission introduces a positively charged hole into 
an insulating lattice, and in the case of Ca2Cu02Cl2, one 
with a highly ionic CaCl rocksalt layer. Although lattice 
polarons have been discussed in the cuprates for some 
time (for instance, see Refs. Il^ and [23l). we believe this 
work represents direct evidence from ARPES of lattice 
polaron formation in the undoped cuprates. 

In simple polaron models, the separation, A£^, between 
the lowest lying state (i.e. the actual QP) and the cen- 
troid of the main peak provides an estimate of the average 
number of phonons, (n), dressing the polaron (llj . For an 
insulator, we take the longitudinal optical (LO) phonon 
branch here to be the characteristic frequency, ^ 40 
meV, yielding (n) = 0.42/0.04 - 10 for the LHB. For 
(n) = 10, the Poisson distribution is effectively indistin- 
guishable from a Gaussian, consistent with our fits. Here 
we use the value of 40 meV, as this is in the range of ener- 
gies for LO phonons established by both optical reflectiv- 
ity and inelastic neutron scattering in La2Cu04 [2^. 12^. 
In principle, choosing loq can significantly change the ab- 
solute value of the coupling constants and (n), but should 
not change our qualitative arguments, at least while re- 
maining within the strong coupling regime. This sur- 
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prisingly large value of (n) puts us well within the small 
polaron limit where the QP residue Z — > and the pho- 
tohole is completely localized. The ARPES spectra then 
consists of a superposition of narrow bands separated by 
ujQ, which in the solid state would be blurred into a single 
Gaussian envelope, as we have observed. Moreover, re- 
cent detailed numerical calculations of the single-particle 
spectral function by Rosch et al. which utilize a, t — J 
model framework to describe the electronic part of the 
system and a realistic shell model for the phonons, also 
suggest a small polaron-like lineshape and broadening 
which agree with what is observed experimentally [26| . 

To further corroborate this, we measured the tempera- 
ture dependence of the LHB, along with the O 2^^^ band 
and the Ca ip core levels, as shown in Fig. [5] where the 
spectra have been scaled and shifted to the peak position 
for the purpose of comparison. The temperature depen- 
dence at the top of the LHB and the top of the O 2p7r 
band are shown with fits in Figs. [2^ and[2)3, respectively. 
Plhb exhibits a much stronger temperature dependence 
than would be expected from conventional electronic in- 
teractions, which in this case would essentially unobserv- 
able, as also reported in Ref. H^. However, the O 2p7r 
peak also exhibits substantial temperature broadening, 
again suggesting that FCB is not unique to the LHB. 

Data from the Ca ip core level taken with He H pho- 
tons are also shown in Fig. [^t after the subtraction of a 
Shirley-like background. Upon closer inspection, the Ca 
2>p core levels are comprised of two spin-orbit split dou- 
blets [1^, resulting in the small shoulder at ^ -17.8 eV. 
We ascribe the second doublet to a surface core level shift 
(SCLS) of Ca^+ in the ionic CaCl cleavage plane arising 
from the altered Madelung potential from the surface ter- 
mination. This has been observed in many materials, in- 
cluding other cuprates such as YBa2Cu307_5 (2^, and 
does not necessarily affect the relevant electronic states 
in the underlying Cu02 plane. Through a simple analysis 
using the known structure and taking the escape depth 
as 10 A ii, we estimate the surface intensity fraction to 
be Sf = 31%, roughly consistent with experiment. Even 
without fitting, the raw data exhibit observable tempera- 
ture broadening, as evidenced by the dip at -19.5 eV (ar- 
row) becoming more pronounced at lower temperatures. 
Fits to a pair of Gaussian doublets are shown, where 
the parameters Sf (33%), spin-orbit splitting (1.19 eV), 
and the SCLS (0.61 eV) were kept fixed with tempera- 
ture. The better agreement of Gaussians for the Ca ip 
core level, as opposed to Lorentzian or Doniach-Sunjic 
lineshapes, again suggests that lattice effects may be the 
dominant broadening mechanism [s^l- 

Before proceeding further, we will rule out spurious 
extrinsic effects which could contribute to the observed 
broadening. Surface roughness was determined to be 
unimportant from well-defined laser reflections (±1°) 
and highly k-dependent valence band spectra. The 
strong temperature and momentum dependence (later 
discussed) of F also rule out the possibility of impurity 
states or electrostatic charging as significant factors. In 
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FIG. 2: Temperature dependence of the LHB at (7r/2, 7r/2) (a) 
and O 2p,r state at (tt, tt) (b). (c) Temperature dependence of 
the Ca 3p core level, with fits to a pair of Gaussian doublets. 
B and S denote the individual contributions of the bulk and 
surface, respectively. In (a)-(c), raw data points are shown 
as symbols, while Gaussian fits are shown as solid lines, (d) 
Temperature dependence of the FWHM of the LHB (circles), 
O 2p^ (squares), and Ca 3p core level (diamonds). 



particular, peaks become sharper at lower temperatures, 
opposite to charging effects which should increase as sam- 
ples become more insulating. We attempted to minimize 
the effects of charging by freshly cleaving the bottom 
side of the sample before mounting in order to minimize 
contact resistance, as well as using the thinnest possi- 
ble samples. Still, charging was typically observed below 
170 K and manifested itself as a higher binding energy 
shift and a broadening of the spectra upon varying the 
incident photon flux. We cannot rule out that some inho- 
mogeneous charging may have contributed to some of the 
sample-to-sample variations in F, represented by the er- 
ror bars in Fig. [JJl. However, we observed similar dT/dT 
on multiple individual cleaves. A breakdown of the sud- 
den approximation or extrinsic losses were ruled out by 
varying the kinetic energy of the photoelectrons between 
9 and 27 eV, while Flhb remained invariant. 

A summary of the temperature dependence is shown 
in Fig. Hi. Between 200 < T < 400 K, dT/dT was es- 
timated to be 1.0 ± 0.3 meV / K for the LHB, 0.6 ± 0.4 
meV / K for the O 2p„, and 0.7 ± 0.2 meV / K for 
the Ca 3p core level. Core level measurements of alkali 
halides such as KCl, where polaron formation is also ob- 
served, give comparable values of dT/dT 0.8 meV / K) 
. The temperature dependence of the FCB can be ex- 
plained as arising from the additional thermal population 
of phonons in a Bose-Einstein distribution. The absolute 
values of F for the three states are also found to be rather 
different. The fact that one observes very different val- 
ues of F and temperature dependences for the LHB, O 



4 



2p7r, and Ca 3p states raises a number of important is- 
sues. The first is that while the lattice contributes to the 
observed broadening of all three states, their coupling 
strengths to various phonon branches could potentially 
be very different depending on the orbital character and 
symmetry of the electronic wavef unctions. Therefore, the 
LHB, O 2p^, and Ca 3p states are likely coupled to differ- 
ent phonon branches with different symmetries, frequen- 
cies, and coupling strengths. This is corroborated by 
calculations in Ref. [2^ which find rather different spec- 
tral distributions for the coupling of phonons to singlet 
states as opposed to non-bonding O 2p orbitals. 

The second issue to be considered is the presence of ad- 
ditional decay channels, in addition to electron-lattice, 
which may also be relevant. For the O and Ca 
3p states which sit deeper in binding energy ( ~ -2 eV 
and -19 eV, respectively), there are interband electron- 
hole decay channels which will contribute to the lifetime 
broadening. Still, the pronounced temperature depen- 
dence observed in the O 2j>^ and Ca ip linewidths in Fig. 
[5] and good agreements to a Gaussian lineshape clearly 
suggest strong electron-phonon contributions. For the 
case of the LHB, the creation of such electron- hole pairs 
is suppressed by the presence of the charge-transfer gap, 
and will therefore not contribute to the linewidth. Still, 
strong electronic / magnetic correlation effects could also 
potentially play an important part in contributing to the 
observed linewidth, in concert with electron-phonon in- 
teractions. As shown in Fig. [21 the LHB clearly ex- 
hibits the strongest temperature broadening of all three 
states, and therefore some enhanced broadening could 
arise due to electronic effects, in combination with lat- 
tice polaron formation. Indeed, theoretical calculations 
indicate that in such strongly interacting systems, the 
electron-electron and electron-phonon interactions can- 
not be considered independently, but can instead feed- 
back upon and enhance one another [sil . [33 | . 

While we have demonstrated the importance of lattice 
effects to the lineshape, strong electron correlations ob- 
viously remain crucial for describing many of the global 
features such as the dispersion of the LHB peak which 
still tracks the extended t — J model calculations 0] , and 
the transfer of spectral weight to higher energies [ssj . 
While the apparently dispersive nature of the LHB might 
appear contradictory to a polaronic scenario, recent the- 
oretical works have resolved this apparent paradox. Cal- 
culations which have incorporated electron-phonon inter- 
actions into the t~ J framework show that while the true 
QP is self-trapped and has vanishing intensity, much of 
the remaining spectral weight forms an incoherent peak 
which tracks the underlying dispersion of the bare t — J 
model [s^l, much like experiment. Other recent work has 
shown in general for polaronic systems that the centroid 
of incoherent weight tracks the dispersion predicted in 
the absence of electron-phonon interactions [131 . This 
also suggests that the t — J model still provides a good 
approximation for the dispersion of the single hole in the 
"frozen" lattice. 
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FIG. 3: (a) Data from the LHB taken at various k points, to- 
gether with Gaussian fits (red), k-dependent dispersion of the 
LHB, £(k) (blue circles), along the (0, 0) — (vr, tt) (b) and (0,7r)- 
(7r,0) (c) directions. Red crosses represent e(k) + 1.8 rLHB(k). 
Data in (a) (c) were taken at 200 K. (d) Comparison of the 
low-energy spectral weight from Cai.95Nao.o5Cu02Cl2 with 
rLHB(k)^^ as a function of angle, as defined in the inset. 



Our experiments have also revealed an intriguing mo- 
mentum anisotropy of Flhb which may have profound 
implications for the emergence of the first doped hole 
states. In particular, our analysis reveals that as a func- 
tion of momentum, the width of the LHB peak, rLHB(k), 
increases directly proportional to the binding energy of 
the peak maximum, £(k), such that e(k) — a Flhb (k). 
Empirically, we have found that a = 1.8, and at 1.8 F 
away from the peak, the intensity has decayed to ^ 10"'* 
of the peak value. In Fig. [5^,, we show data from 
Ca2Cu02Cl2 together with Gaussian fits at various lo- 
cations in k-space, along with crosses to denote the po- 
sition of e(k) -I- 1.8FLHB(k). The momentum depen- 
dence of these features along both the (0, 0) — (tt, tt) and 
(0, tt) — (tt, 0) directions is summarized in Figs. [3)d and 
12};. While the peak position changes considerably, the on- 
set remains roughly fixed, and thus the width appears to 
change to accommodate the dispersion of the peak. Here 
we should note that the simple Gaussian fit continues 
to work well over an extended range in k-space. This is 
somewhat surprising, given that the original assumptions 
for a Poisson distribution (and hence a Gaussian profile) 
for the spectral function were based on the idealized case 
of a single electronic state in a completely nondispersive 
band. Given that the dispersion of the LHB in the frozen 
lattice is ^ 300 meV, this is considerably larger than the 
characteristic phonon frequencies, and therefore it is not 
clear that a Gaussian lineshape would remain a realistic 
approximation. Nevertheless, from an empirical stand- 
point, a Gaussian profile appears to describe the low- 
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energy portion of the spectral lineshape of not only the 
LHB, but also the O 2pjr band quite well. 

In a simple picture where F is related to the electron- 
phonon coupling strength, the photohole would effec- 
tively be coupled more strongly to the lattice at k = 
(tt, 0), the d-wave antinode, than at k = (7r/2,7r/2), the 
d-wave node. Recent results from Ca2-2;Naa:Cu02Cl2 
have demonstrated that spectral weight first emerges 
near (7r/2,7r/2) upon hole doping, while only faint in- 
tensity is visible near {7r,0) [35l. [s^. and similar re- 
sults were reported for other lightly hole doped cuprates 
such as La2-a;Sra;Cu04 and Bi2Sr2CaCu208+5 [33, HI]. 
The momentum dependence of the inverse peak width 
in undoped Ca2Cu02Cl2, rLHB(k)^^, exhibits a qualita- 
tive similarity to the low energy spectral weight in the 
lightly doped materials. This is shown in Fig. where 
the spectral weight along the ostensible Fermi surface 
for Cai.95Nao.o5Cu02Cl2 is integrated in a narrow win- 
dow around Ep (Ep ±10 meV) 36] and compared with 
rLHB(k)~^ along the antiferromagnetic zone boundary 
(dashed red) from the undoped system. The correspon- 
dence between these two quantities may suggest that the 
lack of well-defined QPs near (tt, 0) in the lightly doped 
compounds may be related to the apparently anisotropic 
coupling in the parent compound. 

In conclusion, we have established through a compar- 
ison of the LHB with the O 2pn and Ca 3p states that 



lattice interactions constitute a significant contribution 
to the observed FCB and the photohole exhibits behav- 
ior similar to a small polaron. We also find a com- 
pelling correspondence between the anisotropy of the 
photoholc-lattice coupling and the momentum depen- 
dence of the low-lying spectral weight in the lightly hole 
doped cuprates. While strong correlation effects cer- 
tainly dictate many important aspects, such as the ex- 
istence of the Mott gap and the dispersion of spectral 
weight, our work demonstrates that the inclusion of lat- 
tice interactions should be important to fully explain the 
dynamics of the single hole in the antiferromagnetic insu- 
lator, and we view these results as complementary, and 
not contradictory to results from purely electronic (i.e. 
i — J or Hubbard) models. This work emphasizes that 
realistic calculations of many strongly correlated systems, 
including poorly conducting transition metal oxides and 
Mott insulators, should include the effects of lattice in- 
teractions together with electronic correlations. 
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